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HIGHLIGHTS 


►  Time-dependent  mechanical  behavior  of  reinforced  PFSA  characterized  via  experiments. 

►  Mechanical  response  of  a  composite  membrane  is  investigated  numerically. 

►  Effect  of  humidity  cycling,  thickness,  and  clamping  pressure  is  studied. 

►  Anode  side  reinforced  configuration  develops  lowest  membrane  stresses. 
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The  mechanical  response  of  a  composite  fuel  cell  membrane,  made  from  layers  of  reinforced  and 
unreinforced  PFSA  material,  is  investigated  via  both  experimental  and  numerical  means.  First,  the  time- 
dependent  mechanical  properties  for  the  reinforced  layers  are  measured  for  a  range  of  environmental 
and  loading  conditions.  A  three-network,  viscoelastic-plastic  constitutive  model  is  developed  to  char¬ 
acterize  the  mechanical  response  of  this  reinforced  membrane  material.  This  constitutive  model  is  then 
used  in  finite  element  simulations  of  a  fuel  cell  unit  (consisting  of  composite  membrane,  electrodes,  gas 
diffusion  layer  and  bipolar  plates)  where  the  effect  of  relative  humidity  (RH)  cycling  on  the  stress 
response  of  the  composite  membrane  is  investigated.  Using  numerical  simulations,  various  layering 
configurations  for  the  composite  membrane  and  different  load  cases  are  studied.  The  investigation 
provides  insight  into  the  stress  response  of  the  membrane  and  suggests  possible  configurations  that  may 
improve  the  effective  membrane  life. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  goal  of  increasing  the  durability  and  lifetime  of  polymer 
membrane  fuel  cells  has  motivated  active  research  and  develop¬ 
ment  into  altering  the  thermal  and  mechanical  properties  of  per- 
fluorosulphonic  acid  (PFSA)  membranes  [1-5].  Specifically,  to 
increase  the  mechanical  strength  of  PFSA  membranes,  various 
reinforcing  agents  have  been  used  such  as  polyvinylidene  fluoride, 
expanded  polytetrafluoroethylene  (e-PTFE)  fibrils  and  other  poly¬ 
mers  [6-13]. 

Several  studies  have  shown  an  increase  in  the  durability  of 
membranes  with  the  addition  of  e-PTFE  reinforcement  [14-18]. 
Cleghorn  et  al.  [14]  conducted  single-cell  life-tests  with 
a  commercial  e-PTFE  reinforced  GORE-SELECT®* 1  series  membrane 
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and  reported  a  low  rate  of  performance  degradation  under 
continuous  H2/air  operating  conditions.  Tang  et  al.  [15]  carried  out 
a  study  comparing  Nation®2  membrane  and  e-PTFE  reinforced  PFSA 
membranes,  and  showed  that  the  swelling  and  shrinkage  stresses 
caused  by  water  uptake  are  lower  in  the  reinforced  membrane  and 
concluded  that  this  was  instrumental  in  their  improved  durability. 

Viscoelastic  models  have  been  used  to  characterize  the  time- 
dependent  mechanical  properties  of  reinforced  GORE-SELECT®-57 
membranes  over  a  wide  range  of  temperatures  and  humidities 
[16,17].  These  studies,  along  with  another  [18]  that  considered 
fatigue  and  creep  effects,  showed  that  the  GORE-SELECT® 
membranes  have  considerably  longer  lifetimes  due  to  the  e-PTFE 
reinforcement. 

Typically,  membrane  failure  in  fuel  cells  is  defined  as  a  critical 
volumetric  rate  of  gas  crossover  from  the  anode  side  to  the  cathode 
side  of  the  membrane.  High  compressive  and  residual  tensile 
stresses  developed  in  the  membrane  under  hygrothermal  cycling 
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during  fuel  cell  operation  are  thought  to  be  responsible  for  causing 
fatigue  loading  [19-25].  This,  in  turn  is  associated  with  propagation 
of  cracks  and  pinholes  across  the  membrane-electrode  assembly 
(MEA),  which  lead  to  failure  [20,26].  Chemical  degradation  also 
plays  a  prominent  role  in  the  failure  process.  In  order  to  investigate 
mechanical-only  failures,  accelerated  mechanical  degradation  test 
protocols  have  been  developed.  These  protocols  typically  involve 
cycling  the  membrane  between  high  and  low  RH  conditions  at 
elevated  temperatures,  in  an  open-circuit  to  eliminate  electro¬ 
chemical  effects,  while  monitoring  the  leakage  of  gas  across  the 
membrane  [27-30].  Gittleman  et  al.  [26]  showed  experimentally 
that  after  several  thousand  such  RH  cycles,  unreinforced  PFSA 
membranes  developed  through-cracks  in  the  membrane  thickness 
direction.  However,  in  reinforced  membrane  samples,  the  cracks 
were  arrested  by  the  reinforcement  layer.  Liu  et  al.  [31]  reported 
that  the  reinforcement  of  PFSA  membrane  with  PTFE  improves  the 
dimensional  stability  and  reduces  the  in-plane  swelling  stress  in 
the  membrane.  The  current  study  aims  to  increase  our  under¬ 
standing  of  the  mechanical  behavior  of  the  e-PTFE  reinforced  PFSA 
membrane  material  and  the  failure  processes  for  a  composite 
membrane  consisting  of  two  or  more  material  layers  of  unrein¬ 
forced  PFSA  (Nafion®211)  and  e-PTFE  reinforced  PFSA. 

In  previous  work  [32-34],  we  characterized  the  time,  tempera¬ 
ture  and  hydration  dependent  properties  of  Nation® 211  membrane 
and  simulated  its  mechanical  response  under  hygrothermal  cycling 
in  a  fuel  cell  unit.  The  simulation  results  suggested  that  large 
(9  tolO  MPa)  residual  tensile  stresses  develop  upon  dehydration, 
which  can  lead  to  damage  in  the  membrane.  In  this  work,  we  explore 
the  effect  of  the  e-PTFE  reinforcement  on  the  membrane’s 
mechanical  response  by  simulating  (via  finite  element  simulations), 
a  composite  membrane  (consisting  of  layers  of  unreinforced  and 
reinforced  membrane  material)  under  hygrothermal  conditions  in 
a  fuel  cell  unit  assembly.  We  use  the  time-dependent  properties  of 
Nation®  membrane,  presented  in  earlier  work  [32]  along  with  the 
properties  for  the  reinforced  membrane  material  established  in  this 
work. 

The  description  of  the  test  setup  and  the  experimental  results 
pertaining  to  the  e-PTFE  reinforcement  material  are  provided  in 
Sections  2.1  and  2.2,  respectively.  A  three-network  viscoelastic- 
plastic  constitutive  model  used  to  characterize  the  mechanical 
response  is  discussed  in  Section  2.3.  Using  the  results  from 
Section  2.3,  the  constitutive  model  is  then  incorporated  in 
numerical  analysis  of  a  composite  membrane  in  a  fuel  cell  oper¬ 
ating  under  hygrothermal  loading  in  Section  3.  Finally,  the  finite 
element  simulation  results  for  various  configurations  of  the 
composite  membrane  are  presented  in  Section  4.  The  in-plane 
stresses  in  the  membrane,  which  have  been  shown  to  be  related 
to  crack  initiation  [24,35]  and  further  degradation  of  the 
membrane,  are  reported  for  selected  conditions.  The  effects  of 
symmetric  RH  cycling  (equal  application  from  the  cathode  and 
the  anode  sides),  asymmetric  RH  cycling,  membrane  thickness 
and  clamping  pressure  are  discussed  for  various  layering  config¬ 
urations  of  the  composite  membrane. 

2.  Mechanical  properties  of  the  reinforced  membrane 
material 

2.1.  Test  setup 

We  evaluated  the  time-dependent  mechanical  properties  for  an 
e-PTFE  reinforced  PFSA,  manufactured  specifically  for  this  study 
and  provided  by  W.  L.  Gore  &  Associates  Inc.  The  membranes  are 
manufactured  in  sheets  of  7  pm  nominal  thickness,  and  the  rein¬ 
forcement  is  uniformly  distributed  throughout  the  thickness  of  the 
membrane.  The  reinforced  membrane  material  has  two  distinct 


in-plane  directions  designated  “machine”  and  “transverse”. 
Therefore,  we  cut  100  mm  by  30  mm  rectangular  test  specimens 
from  the  membrane  sheet  in  both  these  directions. 

The  test  protocol  is  similar  to  that  used  in  our  previous  work 
[33,36],  but  repeated  here  for  clarity.  We  conducted  the  tests  using 
the  test  setup  shown  in  Fig.  1  consisting  of  an  MTS  Alliance™  RT/5 
material  testing  system  fitted  with  an  ESPEC  custom-designed 
environmental  chamber.  Uniaxial  tensile  and  stress-relaxation 
tests  were  conducted  under  controlled  environmental  conditions 
and  at  selected  load  parameters  such  as  displacement  rates  and 
hold  strains. 

Each  specimen  was  aligned  with  the  MTS  machine  axis  and 
clamped  in  a  pair  of  vise-action  grips  to  provide  a  nominal  gauge 
length  of  50  mm  as  determined  by  the  grip  separation  (Fig.  1).  To 
achieve  the  desired  environmental  conditions  in  the  chamber,  the 
temperature  was  first  set  to  the  desired  value  and  allowed  to 
stabilize  and  then  the  humidity  was  slowly  increased  (or  decreased) 
to  the  desired  RH  value  with  the  specimen  slack.  Both  the 
temperature  and  humidity  were  kept  at  the  desired  conditions  for 


Fig.  1.  (A)  MTS  alliance  RT/5  material  testing  system  with  custom-designed  environ¬ 
mental  chamber  from  ESPEC  Inc.,  and  (B)  inside  of  the  chamber  showing  the  specimen 
held  with  a  pair  of  vise-action  grips. 
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at  least  a  half  hour  to  ensure  equilibrium  conditions  before  applying 
tension.  During  this  process,  the  length  of  the  specimen  changes 
due  to  the  thermal  and  swelling  deformations  of  the  membrane. 
Therefore,  before  applying  a  force,  the  crosshead  was  manually 
adjusted  until  the  initial  force  applied  to  the  specimen  was  brought 
to  a  small,  finite  tensile  value  ( ~  0.01  N),  eliminating  the  initial  slack 
caused  by  thermal  and  swelling  expansions.  To  determine  the 
subsequent  strain,  we  took  the  undeformed  length  of  the  specimen 
as  the  original  length  at  ambient  conditions,  plus  the  total 
displacement  of  the  crosshead  corresponding  to  the  change  in 
length  caused  by  the  change  in  environmental  conditions. 

We  tested  three  specimens  each  in  the  machine  and  transverse 
directions  at  each  temperature,  humidity  and  tensile  loading  rate 
or  relaxation  holding  strain  combination.  For  each  specimen,  the 
pre-test  thickness  and  width  were  measured  with  a  micrometer 
and  a  caliper,  respectively,  at  three  locations  along  the  sample 
before  testing.  The  averages  of  these  three  measurements  were 
used  as  the  nominal  dimensions  of  the  sample  under  ambient 
conditions. 

2.2.  Experimental  results 

2.2.1.  Uniaxial  tensile  test 

From  the  first  several  sets  of  tests,  we  found  the  mechanical 
response  of  the  reinforced  membrane  to  be  almost  independent  of 
the  humidity  (Fig.  2A).  Therefore,  the  remaining  tests  were  con¬ 
ducted  at  four  temperatures  (25,  45,  65  and  80  °C)  and  a  relative 
humidity  of  30%.  Two  displacement  rates  of  10  mm  min-1  and 
250  mm  min-1  were  selected,  for  each  of  the  two  loading  directions 
(machine  and  transverse),  giving  16  distinct  conditions.  Three 
specimens  were  tested  at  each  condition  and  the  average  response 
is  reported.  The  acquired  test  data,  in  the  form  of  load  as  a  function 
of  displacement,  are  converted  into  true-stress  as  a  function  of 
true-strain  and  presented  in  Fig.  2A— C  for  selected  conditions. 

We  make  the  following  observations  from  the  test  results.  An 
increase  in  load  rate  results  in  an  increase  in  the  elastic  modulus 
and  the  proportional  limit  stress  (Fig.  2B).  There  is  also  an 
observable  increase  in  the  hardening  modulus  with  an  increase  in 
load  rate  (this  behavior  was  not  seen  for  unreinforced  Nation® 
membranes  [32,33]).  The  elastic  modulus  and  the  proportional 
limit  stress  were  found  to  decrease  with  an  increase  in  temperature 
(Fig.  2C).  However,  the  rate  of  change  is  reduced  as  the  temperature 
is  increased  (80  °C).  The  elastic  modulus,  proportional  limit  stress 
and  the  hardening  modulus  (post-proportional  limit  tangent 
modulus)  for  various  environmental  and  loading  conditions  are 
summarized  in  Fig.  3. 

2.2.2.  Stress-relaxation  text 

As  observed  for  the  tensile  test,  there  was  negligible  effect  of 
humidity  on  the  relaxation  response  (Fig.  4A).  Therefore  we  con¬ 
ducted  the  stress-relaxation  tests  at  the  temperatures  of  25, 45,  65 
and  80  °C  and  a  humidity  of  30%RH.  An  initial  displacement  rate  of 
10  mm  min-1  was  used  to  reach  the  desired  strain.  We  recorded  the 
stress  evolution  of  the  membrane  for  three  holding  strain  levels 
(0.05,  0.1  and  0.2)  for  each  of  the  two  directions  of  loading,  giving 
24  test  conditions. 

The  results  show  that  when  the  samples  are  held  at  a  higher 
constant  strain,  the  stress  equilibrates  at  higher  magnitude  (Fig.  4B) 
and  that  the  overall  stress  levels  and  the  equilibrium  stress 
decrease  with  an  increase  in  temperature  (Fig.  4C). 

2.3.  Constitutive  model  for  the  reinforced  membrane  material 

Here,  we  present  the  constitutive  model  used  to  capture  the 
time-dependent  mechanical  response  of  the  reinforced  membrane. 


True  Strain 


Fig.  2.  Selected  experimental  data  showing  true-stress  vs.  true-strain  results  for 
uniaxial  tensile  test  measurements  for  various  (A)  hydration  levels,  (B)  load  rates,  and 
(C)  temperatures. 
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The  constitutive  material  parameters  are  all  functions  of  tempera¬ 
ture  and  strain  rate  (wherever  applicable)  as  determined  from  the 
experiments. 

We  assume  that  the  total  strain  tensor,  ey,  is  given  as  the  sum 

£ij  =  el  +  ejj  +  ef-,  (ij  =  1,2,3),  (1) 

where  ejj,  ej-  and  ej-1  are  the  swelling,  thermal  and  mechanical  strain 
contributions,  respectively.  The  swelling  and  the  thermal  strain 
formulations  have  been  discussed  in  our  previous  work  [32],  but 
are  summarized  in  Appendix  A.l  and  A.2  for  completeness.  The 
mechanical  strains  (e^1)  are  defined  in  the  context  of  a  three- 
network  viscoelastic-plastic  model,  Fig.  5,  where  one  network 
corresponds  to  the  time-independent,  elastic-plastic  response  and 
the  other  two  networks  correspond  to  the  time-dependent,  elastic- 
viscous  response.  By  definition,  the  total  mechanical  strain,  ej-1,  is 
equal  in  all  three  networks 


25  45  65  80 

Temperature  (°C) 


Temperature  (°C) 


25  45  65  80 

Temperature  (°C) 


Fig.  3.  (A)  elastic  modulus,  (B)  proportional  limit  stress,  and  (C)  hardening  modulus  for 
the  reinforced  membrane  at  various  temperatures  for  two  load  rates  (10  and 
250  mm  min-1)  and  two  directions  of  loading  (M-Machine  and  T-Transverse). 


A 


Time  [s] 


Fig.  4.  Sample  experimental  data  showing  true-stress  vs.  time  evolution  for  various 
(A)  hydration  levels  (at  two  hold  strains  of  0.05  and  0.1),  (B)  holding  strains  (at 
a  temperature  of  25  °C),  and  (C)  temperatures  (at  a  hold  strain  of  0.2). 
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EP  :  Elastic-Plastic 
EV!  :  Elastic-Viscous  1 
EV2 :  Elastic-Viscous  2 

EV2  EVt  EP 


N - H  M - H 

Time  Time 

Dependent  /^dependent 


Fig.  5.  One-dimensional  idealization  of  the  three-network  viscoelastic-plastic  consti¬ 
tutive  model. 


M 


_  „EP  _  tvi  _ 


EV, 


EV2 


(2) 


Fig.  6.  True-stress  vs.  true-strain  uniaxial  tensile  response  of  the  reinforced  membrane 
material  at  various  conditions  obtained  from  constitutive  model  predictions  compared 
with  the  experiments. 


where  subscript  el  refers  to  the  elastic  portion  and  vi  refers  to  the 
viscous  portion  of  the  strain.  Assuming  a  Neo— Hookean  hypere¬ 
lastic  formulation  for  the  constitutive  response  of  the  spring  in  this 
network  [37],  the  Cauchy  stress,  a^2t  is  given  by 


where  the  superscripts  EP,  EVi  and  EV2  refer  to  the  elastic-plastic, 
the  first  elastic-viscous  network  and  the  second  elastic-viscous 
network  of  the  model,  respectively  (Fig.  5).  The  strain  formula¬ 
tions  for  the  elastic-plastic  (EP)  network  and  the  first  elastic- 
viscous  (EVi)  network  are  discussed  in  our  previous  work  [32] 
and  summarized  in  Appendix  A.3  and  A.4,  respectively,  while  the 
following  discussion  describes  the  strains  in  the  elements  of  the 
second  elastic- viscous  (EV2)  network. 

The  total  strain  in  the  second  elastic- viscous  (EV2)  network  is 
the  sum  of  the  strain  in  the  spring  and  the  dashpot  element  given 
by 


(3) 


Table  1 

(A)  elastic-plastic  and  (B)  elastic-viscous,  constitutive  model  parameters  for  the 
e-PTFE  reinforced  PFSA  membrane  material. 


A 

Temp 

(°C) 

Kep 

(MPa) 

ay  H 

(MPa)  (MPa) 

/<eVi  (MPa) 

I<ev2 

(MPa) 

10  mm  min-1 

250  mm  min  1 

25 

200 

5 

150 

539 

622 

8.3 

45 

120 

4 

100 

130 

274 

4.8 

65 

100 

4 

72 

111 

211 

3.9 

80 

86 

3 

60 

100 

204 

3.2 

B 

Temp 

In  (A) 

n 

B  m 

(°C) 

10  mm  min  1 

250  mm  min 

1  10  mm  min 

1  250  mm  min  1 

25 

-23.03 

-23.03 

8 

7 

50  6.4 

45 

-20.72 

-23.03 

8 

7 

50  5.5 

65 

-17.32 

-19.11 

8 

6.8 

50  3.3 

80 

-18.42 

-20.72 

8 

6.5 

50  2.6 

af2  =  ^  dev  [by]  +  k(J  -  1  )<5„ ,  (4) 

where  /i  =  I<ev2/ 2(1  +  v)  and  k  =  I<Ey 2/3(l  -  2v)  are  the  shear  and 
bulk  modulus  respectively,  with  KEy2  being  the  elastic  modulus.  The 
tensor  by  is  the  distortional  left-Cauchy  Green  tensor  given  by 
3FikFjk  and  the  Jacobian,  J ,  is  given  by  det(F),  where  F  is  the  defor¬ 
mation  gradient. 

The  dashpot  in  this  network  follows  Hill’s  criteria  [38]  to 
model  anisotropic  viscous  behavior.  The  generalized  three- 


Time  [s] 


Fig.  7.  True-stress  evolution  as  a  function  of  time  for  the  uniaxial  stress-relaxation 
test  of  the  reinforced  membrane  material  at  various  conditions  obtained  from 
constitutive  model  predictions  compared  with  the  experiments. 
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Finite  Element  Model 


Fig.  8.  Schematic  of  the  2D  generalization  of  a  single  fuel  cell  unit  used  to  conduct  composite  membrane  simulations  (not  to  scale). 


[A]  Unreinforced  PFSA 


[B]  Core  reinforced 


u-PFSA 

|6.25|im 

25|im 

u-PFSA 

25nm 

r-PFSA 

J12.5nm 

u-PFSA 

t6.25nm 

[C]  Outer  reinforced 

[D]  Anode  side  reinforced 

'888  r-PFSA 

t6.25|im 

r-PFSA 

Jl2.5nm 

25nm 

u-PFSA 

Jl2.5nm  25 nm 

t 

■  r-PFSA 

t6.25um 

u-PFSA 

|l2.5|im 

Fig.  9.  Schematic  configurations  of  the  composite  membrane  with  various  arrangements  of  the  reinforcement  layer.  (A)  unreinforced  PFSA  (reference),  (B)  core-reinforced,  (C) 
outer-reinforced,  and  (D)  anode  side  reinforced  configuration.  u-PFSA  and  r-PFSA  refer  to  unreinforced  and  reinforced  membrane  material,  respectively. 
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Table  2 

Material  properties  for  various  fuel  cell  components  used  in  the  finite  element  model. 


Behavior 

Properties 

Swelling  behavior 

Thermal  expansion 

Reinforced  membrane 

Viscoelastic- 

plastic 

Section  2.3 

Transverse  isotropic  (Ct-of-pUne  =  10C-Piane) 

1.2  •  10-8  K-1 

Unreinforced  membrane 

Viscoelastic- 

plastic 

Viscoelastic- 

[32] 

Isotropic  (cQut_0f_pjane  —  f'jn  p|Jne ) 

1.2-10  8  I<  1 

Catalyst 

[47] 

Assumed  not  to  swell 

Neglected 

plastic 

GDL 

Transverse 

En  =  1500  MPa 

Assumed  not  to  swell 

Neglected 

isotropic,  linear- 
elastic 

E22  =  9  MPa  [48] 

Bipolar  plate 

Linear-elastic 

E  =  10  GPa,  v  =  0.25  [25] 

Assumed  not  to  swell 

Neglected 

dimensional  form  of  the  flow  rate  for  the  dashpot  in  this  network 
is  given  by 


where  m  is  the  flow  exponent  and  B  is  the  flow  resistance  in  this 
network.  The  Hill  stress,  (jHl11  is  given  by 


constructed  a  model  for  a  representative  repeating  element 
adapted  from  our  previous  work  [25,32,35,40]  consisting  of  a  half  of 
a  land  and  a  half  of  a  groove  with  the  dimensions  given  in  Fig.  8. 
Generalized  plane  strain  conditions  are  assumed,  corresponding  to 
a  uniform  strain  in  the  z-direction  (out-of-plane  in  the  figure).  The 
right  edge  and  bottom  have  symmetric  boundary  conditions 
imposed.  To  account  for  possible  displacement  in  the  in-plane  (x) 


—  \J £(@22  ~  cr33)2+G((733  —  d"||)2+H(( 7\\  —  (722)2  +2L(7^3  +  2M(7^1  +  2N(7^2 


(6) 


where  F,  G,  H,  L,  M  and  N  are  the  Hill  parameters  determined  by 
experiment. 

The  parameters  of  the  three-network  viscoelastic-plastic 
constitutive  model  discussed  above  were  obtained  from  the 
experimental  data  for  the  reinforced  membrane  material,  and  are 
given  in  Table  1.  To  verify  that  these  material  parameters  corre¬ 
spond  to  the  mechanical  properties  of  the  membrane,  finite 
element  simulations  of  a  membrane-strip  model  using  the  above 
discussed  constitutive  behavior  were  compared  with  the  experi¬ 
mentally  obtained  stress-strain  and  relaxation  results  (Figs.  6  and 
7).  A  plane  stress  ABAQUS  [39]  FE  model  was  developed  with 
dimensions  conforming  to  the  test  specimen.  The  simulation 
matches  the  experimental  results  fairly  well  and  captures  key 
features  of  the  tensile  response  such  as  elastic  modulus,  onset  of 
non-linearity,  and  hardening  modulus  (Figs.  6  and  7).  In  the  stress- 
relaxation  results,  the  simulation  captures  the  stress  peak,  relaxa¬ 
tion  time  and  the  equilibrium  stress  for  the  observed  time  frame  of 
about  1—1.5  h.  The  trend  suggests  that  relaxation  continues  beyond 
the  test  time,  although  at  a  decreasing  rate.  (Due  to  limitations  of 
the  testing  equipment  and  practical  considerations,  relaxation 
beyond  1.5  h  was  not  measured.) 

3.  Numerical  implementation 

In  this  section,  the  mechanical  response  of  a  composite  PFSA 
membrane  in  a  fuel  cell  subjected  to  hygrothermal  cycling  is 
investigated  via  finite  element  simulations.  We  used  the  commercial 
software  ABAQUS™  [39]  with  8-noded  biquadratic  coupled 
temperature-displacement  generalized  plane  strain  elements 
(CPEG8T).  The  model  consists  of 32,604  elements  and  101,424  nodes 
with  20  elements  along  the  thickness  direction  of  the  composite 
membrane.  Numerical  convergence  was  verified  using  progressively 
finer  meshes  until  no  appreciable  change  was  observed  in  the 
results. 

3.1.  Geometry 

A  typical  fuel  cell  unit  contains  repetitive  grooves  and  lands 
forming  the  gas  channels  in  the  bipolar  plate.  Therefore,  we 
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Fig.  10.  (A)  equilibrium  water  mass  uptake  for  the  unreinforced  and  the  reinforced 
membrane  material  at  80  °C.  (B)  in-plane  length  change  (%)  and  (C)  length  change  after 
normalizing  by  the  final  equilibrium  length,  for  both  the  materials.  The  data  is  based 
on  Griffith  et  al.  [41  ]. 
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Fig.  11.  Transversely  isotropic  swelling  strain  in  the  reinforced  (r-PFSA)  membrane 
material  and  isotropic  swelling  strain  in  the  unreinforced  (u-PFSA)  membrane  material 
as  functions  of  water  volume  fraction. 

direction,  a  continuous  boundary  condition  is  applied,  under  which 
all  the  nodes  on  left  edge  of  the  model  can  move  laterally  in  unison. 
A  uniform  pressure  of  1  MPa  is  applied  on  the  top  edge  of  the 
bipolar  plate  to  simulate  a  typical  spring-loaded  clamping  pressure 


Time  [s] 


Time  [s] 


Fig.  12.  (A)  symmetric  RH  cycling  test  profile  and  (B)  asymmetric  RH  cycling  profile. 
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Fig.  13.  Maximum  tensile  and  maximum  compressive  stresses  (MPa)  in  the  in-plane 
and  the  out-of-plane  directions  for  the  unreinforced  and  the  reinforced  membrane 
materials. 


for  a  fuel  cell.  The  effect  of  altering  the  clamping  pressure  is  also 
investigated  parametrically  in  Section  4. 

The  composite  membrane  is  modeled  as  a  layered  structure 
consisting  of  two  materials  (unreinforced  PFSA  and  reinforced 
PFSA)  arranged  in  four  configurations  (Fig.  9): 

i.  Unreinforced  membrane  (u-PFSA)  material  only  (reference 
case)  (Fig.  9A). 

ii.  Reinforced  (r-PFSA)  membrane  material  as  the  core  sur¬ 
rounded  by  two  layers  of  unreinforced  membrane  material  on 
the  outside  (Fig.  9B). 

iii.  Unreinforced  membrane  material  as  the  core  surrounded  by 
two  layers  of  reinforced  membrane  material  on  the  outside 
(Fig.  9C). 

iv.  Reinforced  membrane  material  on  the  anode  side  and  unre¬ 
inforced  membrane  material  on  the  cathode  side  (Fig.  9D). 

During  manufacturing,  the  various  layers  are  produced  together 
using  a  casting  process.  Therefore,  the  layers  are  connected  in  the 
fuel  cell  model  using  a  ‘tie  constraint’  [39]  in  which  the  interface 
nodes  are  constrained  to  move  together  (i.e.  no  interfacial 
debonding). 
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■  observation  element 
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■  observation  element 
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Fig.  14.  Symmetric  RH  cycle:  In-plane  stress  evolution  for  an  element  each  on  the 
groove  side  in  the  unreinforced  and  the  reinforced  layers  of  the  core-reinforced 
composite  membrane  configuration  for  ten  cycles.  u-PFSA  and  r-PFSA  refer  to  the 
unreinforced  and  the  reinforced  membrane  materials,  respectively. 
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3.2.  Properties 

The  material  properties  used  for  the  various  fuel  cell  compo¬ 
nents  in  the  finite  element  model  are  summarized  in  Table  2.  The 
membrane  is  the  only  component  that  is  assumed  to  swell.  The 
sorption  and  swelling  behavior  used  in  the  model  of  the  membrane 
is  described  in  the  following  section. 

3.2.1.  Sorption  behavior 

The  water  sorption  and  the  consequent  swelling  behavior  of  the 
membrane  used  in  this  work  is  based  on  the  equilibrium  water 
mass  uptake  and  in-plane  swelling  data  measured  by  Griffith  et  al. 
[41].  In  this  experiment,  the  membranes  were  exposed  to  an 
outside  humidity  change  from  dry  (A  =  2)  to  saturated  (A  =  14) 
condition  (at  a  constant  temperature  of  80  °C),  while  the  increase  in 
mass  and  change  in  length  were  measured.  The  data  (Fig.  10A) 
show  that  the  e-PTFE  reinforcement  does  not  significantly  affect 
the  total  water  mass  uptake  [41  ].  Flowever,  dimensional  changes  in 
the  in-plane  direction  are  significantly  different  as  shown  by  the  in¬ 
plane  swelling  measurements  [41]  (Fig.  10B).  The  reinforced 
membrane  swells  considerably  less  in  the  in-plane  direction  than 
the  unreinforced  membrane.  (The  out-of-plane  swelling  can  be 
approximated  by  considering  the  total  water  uptake  of  the  rein¬ 
forced  membrane;  Fig.  10A.)  After  normalizing  by  the  equilibrium 
length,  the  swelling  kinetics  of  the  two  membranes  appears  to  be 


Fig.  16.  Symmetric  RH  cycle:  In-plane  stress  evolution  over  two  cycles  in  the  groove 
side  elements  for  the  unreinforced  (u-PFSA)  and  reinforced  (r-PFSA)  layer  of  the 
composite  membrane.  CR,  OR  and  ASR  refer  to  core-reinforced,  outer-reinforced  and 
anode  side  reinforced  configurations,  respectively. 
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B.  Maximum  in-plane  stress  after  hydration 
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Fig.  15.  Symmetric  RH  cycle:  Snapshots  of  maximum  in-plane  stress  contours  in  the  composite  membrane  after  (A)  clamping,  (B)  hydration,  and  (C)  dehydration  step,  (i)-(iv)  refer 
to  the  four,  layered  configurations  as  shown  in  Fig.  10  u-PFSA  and  r-PFSA  refer  to  the  unreinforced  and  the  reinforced  membrane  materials,  respectively. 
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Time  [s] 

Fig.  17.  Symmetric  RH  cycle:  In-plane  swelling  strain,  e™,  and  the  total  plastic  strain 
magnitude,  ep,  evolution  for  an  element  from  the  groove  side  each  in  the  unreinforced 
(u-PFSA)  and  the  reinforced  (r-PFSA)  layer  of  the  composite  membrane. 


Normalized  distance  from  mid-channel  to  mid-land 


Fig.  18.  Symmetric  RH  cycle:  Distribution  of  the  in-plane  stress  after  dehydration  as 
a  function  of  the  normalized  distance  from  mid-channel  to  mid-land  for  three  values  of 
clamping  pressure:  0.2, 1.0  and  2.0  MPa.  u-PFSA  and  r-PFSA  correspond  to  the  unre¬ 
inforced  and  the  reinforced  layers  of  the  composite  membrane. 


similar  in  nature  as  shown  in  Fig.  10C.  Therefore,  in  the  numerical 
analyses  that  follow,  sorption  kinetics  of  the  reinforced  membrane 
and  unreinforced  membrane  are  assumed  to  be  similar. 

The  experiments  show  that  the  in-plane  swelling  in  the  rein¬ 
forced  membrane  is  about  a  quarter  of  that  of  the  unreinforced 
membrane  [41  ].  If  we  assume  that  the  swelling  of  the  unreinforced 
membrane  is  isotropic,  then  the  out-of-plane  swelling  of  the  rein¬ 
forced  membrane  is  approximately  ten  times  larger  than  its  in¬ 
plane  swelling  as  shown  in  Fig.  11,  to  account  for  the  nearly  equal 
total  water  uptake.  Therefore,  we  assumed  transversely  isotropic 
swelling  behavior,  with  out-of-plane  swelling  equal  to  ten  times 
than  that  of  in-plane  swelling,  for  the  reinforced  material  in  the 
finite  element  simulations. 


3.3.  Load  cases 

A  relative  humidity  (RH)  cycling  test  has  been  developed  by  W.  L. 
Gore  [28]  as  an  accelerated  test  protocol  to  assess  the  mechanical 
durability  of  fuel  cell  membranes.  This  cycle  operates  at  a  constant 
temperature  of  80  °C,  while  the  humidity  is  cycled  between  dry 
(RH  ~  15%)  and  supersaturated  (corresponding  to  wet  gas  dewpoint 
of  90  °C)  vapor  conditions.  The  hygro thermal  cycling  load  cases 
investigated  in  this  work  are  based  on  this  cycle,  referred  to  as  the 
“Gore  cycle”.  In  particular,  two  load  cases  are  considered  in  this 
work: 

3.3.1.  Symmetric  RH  cycle 

The  standard  Gore  cycle,  with  simultaneous  application  of 
equally  humidified  nitrogen  on  both  the  anode  and  the  cathode  side 
gas  channels  is  referred  to  as  ‘Symmetric  RH  Cycle’  in  this  work 


Table  3 

Symmetric  RH  cycling:  Maximum  in-plane  residual  tensile  stress  (in  MPa)  in  the 
unreinforced  (u-PFSA)  and  the  reinforced  (r-PFSA)  membrane  material  for  various 
thicknesses  of  u-PFSA  with  a  constant  reinforcement  core  thickness  (5  pm). 


Thickness  of  the  unreinforced  segment 

1  pm  5  pm  10  pm  20  pm 

u-PFSA  10.38  10.27  10.13  085 

r-PFSA  0.04  0.04  0.03  0.01 


(Fig.  12A).  The  maximum  water  volume  fraction,  0,  corresponding  to 
saturated  vapor  conditions  is  0.32.  Using  equation  (A-2),  this 
corresponds  to  water  content  in  the  membrane  (A)  of  approximately 
14.  The  water  volume  fraction  for  the  dry  state  corresponds  to 
0  =  0.06  (or  A  =  2).  The  loading  is  applied  at  the  GDL  interface  in  the 
channels  for  10  s  corresponding  to  the  wet  state  and  50  s  for  the  dry 
state.  The  transition  from  dry  to  wet  and  vice-versa  is  assumed  to  be 
linear  and  takes  place  in  a  very  short  time  (1  s  in  the  simulation). 


3.3.2.  Asymmetric  RH  cycle 

In  a  working  fuel  cell  drawing  current,  equal  saturation  at  the 
anode  and  the  cathode  is  unlikely.  The  generation  of  water  tends  to 
keep  the  cathode  side  hydrated  and  the  electro-osmotic  drag  tends 
to  dry-out  the  anode  side  with  increasing  current  density  [42]. 
Therefore,  an  ‘Asymmetric  RH  Cycle’  was  investigated  which 
involves  different  RH  loading  at  the  anode  and  the  cathode  sides 
flow  channels  thus  producing  a  water  content  gradient  across  the 
membrane  thickness  (Fig.  12B).  In  this  cycle,  the  cathode  side  flow 
channel  is  kept  hydrated  at  all  times  while  humidity  cycling 
following  the  Gore  cycle  (described  above)  is  applied  at  the  anode 
side  flow  channel  (Fig.  12B). 
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Fig.  19.  Water  and  temperature  profile  in  the  membrane  for  asymmetric  RH  loading. 
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In  both  the  load  cycles,  humid  and  dry  nitrogen  is  alternately 
introduced  into  the  flow  channels.  The  actual  water  profile  at  the 
membrane  is  mediated  by  the  diffusion  through  GDL,  electrode  and 
the  composite  membrane  itself.  The  water  volume  fraction  in  the 
membrane  is  modeled  as  a  third  degree  of  freedom  in  addition  to 
the  two  nodal  displacements  in  the  two-dimensional  generalized 
plane  strain  element  used  for  the  finite  element  analysis.  Therefore, 
a  coupled  sorption-displacement  analysis  is  carried  out  in  which 
the  swelling-induced  stresses  and  the  water  uptake  through 
diffusion  are  coupled  and  solved  simultaneously  [43]. 

4.  Results  and  discussion 

The  hydration  load  on  the  membrane  leads  to  volumetric 
swelling  from  the  water  absorption.  However,  the  swelling  is 
limited  due  to  the  mechanical  constraints  from  the  other  fuel  cell 
components  and  consequently,  gives  rise  to  stresses.  In  our 
previous  work  [25,32,43,44]  we  showed  that  the  in-plane  stresses 
are  the  dominant  components  in  the  unreinforced  membrane 
during  fuel  cell  operation.  This  holds  for  the  composite  membrane 


as  illustrated  by  the  relative  stress  magnitudes  in  the  membrane  for 
a  typical  simulation  (Fig.  13).  Consequently,  only  the  in-plane  stress 
component  will  be  discussed  in  the  following. 

4.1.  Symmetric  loading 

The  typical  in-plane  stress  evolution  in  the  unreinforced  and  the 
reinforced  layers  of  the  core-reinforced  composite  membrane 
configuration,  for  ten  cycles  of  symmetric  RH  cycling  is  shown  in 
Fig.  14.  The  evolution  is  presented  for  a  representative  element  in 
each  of  the  two  layers  of  the  composite  membrane  on  the  groove 
side.  While  both  layers  develop  similar  compressive  stresses  after 
hydration,  the  residual  tensile  stress  in  the  reinforced  membrane 
material  is  significantly  lower  compared  to  the  unreinforced 
material.  This  can  be  attributed  to  the  greater  hydration- 
dependence  of  the  material  properties  for  the  unreinforced 
membrane.  The  peak  stresses  increase  for  the  initial  few  cycles  and 
become  constant  thereafter. 

A  more  detailed  understanding  of  the  stress  evolution  may  be 
obtained  by  studying  the  snapshots  of  the  in-plane  stress  contours 
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Fig.  20.  Asymmetric  RH  cycling:  In-plane  stress  evolution  for  three  elements  (shown  at  the  top)  on  the  groove  side  of  the  composite  membrane  for  (A)  unreinforced  PFSA  only 
(reference  case),  (B)  core-reinforced,  (C)  outer-reinforced,  and  (D)  anode  side  reinforced  configurations. 
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Fig.  21.  Asymmetric  RH  cycling:  Snapshots  of  the  in-plane  stress  contours  in  the 
composite  membrane  after  dehydration,  highlighting  the  residual  tensile  stresses 
developed  in  various  configurations  (A-D)  as  shown  in  Fig.  10.  u-PFSA  and  r-PFSA 
correspond  to  the  unreinforced  and  the  reinforced  layers  of  the  composite  membrane. 


in  the  composite  membrane  for  symmetric  RH  cycling  (Fig.  15).  The 
stress  evolution  over  time  for  three  elements  on  the  groove  side  of 
the  composite  membrane  is  also  shown  (Fig.  16).  After  the  clamping 
pressure  (1  MPa)  is  applied,  small  compressive  stresses  appear 
throughout  the  composite  membrane  (Fig.  15 A).  When  the 
temperature  is  increased  from  the  initial  25  °C  to  80  °C  along  with 
the  increase  in  hydration  from  the  initial  dry  state  (0  =  0.06)  to 
saturated  vapor  (0  =  0.32),  high  compressive  stresses  develop  in 
both  layers  of  the  composite  membrane,  with  the  magnitude  being 
higher  in  the  reinforced  layer  (Fig.  15B).  The  region  under  the  land 
has  higher  compressive  stress  than  the  region  under  the  groove. 
The  system  is  then  held  for  10  s  at  the  maximum  hydration  during 
which  relaxation  of  the  in-plane  stress  occurs  (Fig.  16).  When  the 
system  is  dehydrated  back  to  the  dry  state,  the  membrane  shrinks. 
The  unreinforced  layer  develops  large  tensile  stresses  after  dehy¬ 
dration  while  very  small  residual  tensile  stresses  develop  in  the 
reinforced  layer  of  the  composite  membrane  (Fig.  15C). 

Due  to  the  low  in-plane  swelling  and  relatively  larger  propor¬ 
tional  limit  stress,  the  reinforced  membrane  develops  lower  in¬ 
plane  swelling  strain  and  lower  total  plastic  strain  than  the  unre¬ 
inforced  layer  (Fig.  17).  Therefore,  the  reinforcement  layer  will  most 
likely  reduce  (or  delay)  the  formation  of  pinholes  and  through- 
cracks. 


4.1.1.  Effect  of  thickness 

Simulations  for  configurations  having  different  unreinforced 
outer  layer  thickness  (1, 5, 10  and  20  pm)  with  a  constant  reinforced 
core  thickness  (5  pm)  were  investigated.  The  results  for  symmetric 


RH  cycling  show  that  the  resulting  maximum  residual  tensile  stress 
after  dehydration  decreases  marginally  with  increased  thickness  of 
the  outer  layer,  Table  3. 

4.1.2.  Effect  of  clamping  pressure 

The  effect  of  clamping  pressure  on  the  membrane  stresses  is 
investigated  for  three  different  clamping  loads  (0.1, 1.0  and  2.0  MPa) 
applied  on  top  of  the  bipolar  plate.  The  larger  clamping  pressure 
induces  larger  compressive  stresses  during  hydration  in  the 
membrane  and  consequently,  lower  residual  tensile  stress  after 
dehydration  (Fig.  18).  This  suggests  that,  higher  clamping  pressure 
may  have  a  beneficial  effect  on  membrane  life. 

4.2.  Asymmetric  loading 

Fig.  16  shows  that  the  stress  evolution  in  the  composite 
membrane  is  independent  of  the  relative  location  of  the  reinforce¬ 
ment  layer.  This  is  due  to  the  fact  that  the  stresses  result  from  the 
swelling  strains,  which  are  driven  by  equal  RH  loading  from  both 
sides  of  the  membrane.  The  stresses  in  the  various  layers  would 
depend  on  the  layering  geometry  in  the  presence  of:  (i)  any 
geometric  perturbations  (such  as  cracks  or  debonding)  or  (ii) 
asymmetric  RH  loading  from  the  opposite  sides  of  the  membrane. 
While  the  case  of  geometric  perturbations  is  left  for  a  future  study, 
we  investigate  asymmetric  RH  loading  (Section  3.3)  in  the  following. 

The  resulting  water  and  temperature  profiles  in  the  membrane 
are  as  shown  in  Fig.  19  for  an  asymmetric  RH  cycle  applied  at  the 
gas  channels.  The  stress  results  in  Fig.  20  are  provided  for  three 
elements,  each  located  at  4  pm  from  the  left  boundary.  The  y- 
location  of  the  elements  is  3  pm  from  the  membrane— anode 
interface,  mid-plane,  and  3  pm  from  the  membrane-cathode 
interface.  The  results  suggest  that  the  stresses  on  the  anode  side 
of  the  membrane  are  higher  than  on  the  cathode  side  due  to  the 
larger  variation  in  hydration  on  the  anode  side. 

In  all  the  layering  configurations,  the  reinforced  material 
develops  large  compressive  stresses  throughout  after  hydration. 
However,  no  residual  tensile  stresses  develop  in  the  reinforced 
material  after  dehydration  (Fig.  20B-D).  On  the  other  hand,  in  the 
unreinforced  layer  of  the  core-reinforced  configuration  (Fig.  20B), 
the  maximum  residual  tensile  stress  is  similar  to  the  stress 
magnitude  in  the  reference  configuration  (Fig.  20A).  This  can  be 
explained  by  the  fact  that  the  unreinforced  material  is  on  the 
outside  of  the  composite  membrane  and  is  therefore  exposed  to 
similar  humidity  cycling  as  in  the  reference  case. 

For  the  outer-reinforced  configuration  (Fig.  20C),  the  maximum 
stress  in  the  unreinforced  layer  is  equal  to  the  maximum  stress 
reached  in  the  middle  of  the  membrane  in  the  reference  case.  In  the 
anode-reinforced  configuration,  in  which  the  reinforcement  is  on 
the  anode  side  (Fig.  20D),  the  in-plane  residual  stress  in  the  unre¬ 
inforced  layer  is  similar  to  the  stress  level  on  the  cathode  side  of  the 
membrane  in  the  reference  case.  The  anode  side  reinforced  config¬ 
uration  develops  the  lowest  magnitude  of  residual  in-plane  tensile 
stresses,  which  are  concentrated  on  a  smaller  region  in  the  unrein¬ 
forced  layer  of  the  composite  membrane  (Fig.  21). 

A  summary  of  the  maximum  residual  tensile  stresses  in  the 
unreinforced  and  the  reinforced  layers  of  the  composite  membrane 
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r-PFSA 
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:  RH  cycling 
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Anode  side  reinforced 

T87H 
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Fig.  22.  Maximum  residual  in-plane  stress  (in  MPa)  in  the  unreinforced  (u-PFSA)  and  the  reinforced  (r-PFSA)  layers  of  the  membrane  for  the  investigated  configurations. 
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for  the  various  configurations  is  provided  in  Fig.  22.  The  residual 
tensile  stresses  in  the  unreinforced  material,  in  the  case  of  asym¬ 
metric  cycling,  depend  on  the  relative  arrangement  of  the  two 
membrane  material  layers,  with  the  anode-side  reinforced  config¬ 
uration  leading  to  lowest  maximum  tensile  stress  level  (10%  of  the 
reference  case  maximum  stress). 

5.  Concluding  remarks 


where  6  and  0O  are  the  current  and  reference  temperature  (in  K), 
respectively.  The  term  is  the  anisotropy  ratio  giving  the  dimen¬ 
sional  change  in  the  i-direction  in  the  membrane  relative  to  the 
total  volume  change  of  the  membrane  such  that  +  ?y  +  ?z  =  1.  The 
term  </>p  is  the  polymer  volume  fraction  given  by 


In  this  work,  the  time-dependent  mechanical  properties  of  an  e- 
PTFE  reinforced  experimental  membrane  material  have  been 
investigated  for  a  range  of  environmental  and  loading  conditions 
using  uniaxial  tensile  and  stress-relaxation  tests.  The  constitutive 
response  of  the  reinforced  material  is  characterized  using  a  three- 
network  viscoelastic-plastic  model. 

To  investigate  how  the  reinforcement  would  affect  the 
mechanical  stresses  in  an  operating  fuel  cell  membrane, 
a  composite  membrane  consisting  of  the  unreinforced  PFSA  and  the 
e-PTFE  reinforced  PFSA  material  in  layered  geometries  was  simu¬ 
lated.  The  constitutive  model  for  the  reinforced  membrane  material 
was  incorporated  into  a  finite  element  model  representing  a  single 
fuel  cell  unit  assembly.  The  fuel  cell  FE  analysis  was  conducted  as 
a  coupled  sorption-displacement  problem  [32,43]  enabling  us  to 
model  the  water  content  profile  and  mechanical  stresses  simulta¬ 
neously.  The  stress  evolution  for  various  configurations  of  the 
composite  membrane  as  well  as  the  effects  of  symmetric  RFI 
cycling,  asymmetric  RH  cycling,  membrane  thickness  and  clamping 
pressure  were  investigated. 

The  results  of  the  simulations  show  that  (i)  the  reinforced 
material  in  the  composite  membrane  develops  large  compressive 
stresses  after  hydration  but  significantly  lower  residual  in-plane 
stresses  after  dehydration  when  compared  to  the  unreinforced 
membrane  material;  (ii)  the  relative  arrangement  of  the  rein¬ 
forcement  in  the  composite  membrane  does  not  affect  the 
maximum  in-plane  stress  in  the  unreinforced  layer  during 
symmetric  RH  cycling;  (iii)  the  total  plastic  strain  magnitude  in  the 
reinforced  layer  is  significantly  lower  than  in  the  unreinforced 
layer;  (iv)  increasing  the  thickness  of  the  unreinforced  layer 
decreases  the  residual  in-plane  tensile  stresses  in  both  the  layers  of 
the  composite  membrane;  (v)  an  increase  in  the  clamping  pressure 
leads  to  lower  maximum  residual  in-plane  tensile  stresses;  (vi) 
during  asymmetric  RH  cycling,  the  anode  side  reinforced  configu¬ 
ration  develops  lowest  in-plane  stresses. 

In  conclusion,  these  results  suggest  that  the  addition  of  an  e-PTFE 
reinforcement  material  can  affect  the  stresses,  and  therefore,  the 
mechanical  durability  of  the  membrane.  In  particular,  during  asym¬ 
metric  RH  cycling,  the  residual  in-plane  tensile  stresses,  which  can 
play  a  significant  role  in  mechanical  degradation  and  membrane 
lifetime,  are  lowered  significantly  for  anode  side  reinforced  configu¬ 
ration.  Our  numerical  predictions  therefore,  can  provide  a  basis  for  the 
design  and  optimization  of  composite  PEM  fuel  cell  membranes  with 
the  aim  of  reducing  life-cycle  cost  by  improving  fatigue  resistance. 


EW/pP 

,p  18A  +  EW/pp’ 


(A-2) 


where  A  is  the  number  of  water  molecules  attached  to  each  sul- 
phonic  acid  group,  EW  is  the  equivalent  weight  of  the  ionomer 
(=1100)  and  pp  is  the  density  of  the  dry  ionomer  material 
(=1970  kg  nr3). 


A.2.  Thermal  strains 


Assuming  isotropic  thermal  expansion,  the  thermal  strains  in 
the  membrane  caused  by  a  change  in  temperature  (0  -  0O)  are  given 
by 

el  =  a(6  -  e0)8y,  (A-3) 


where  a  is  the  linear  coefficient  of  thermal  expansion  and  dy  is  the 
Kronecker  delta  given  by 


1;  i=j 

0;  i  =tj  ' 


(A-4) 


A.3.  Strain  in  the  elastic— plastic  (EP)  layer 


Extrapolating  the  spring  concept  in  the  elastic— plastic  layer  to 
represent  a  generalized  three-dimensional  isotropic  Hooke’s  law 
response,  the  elastic  strain  can  be  written  as  a  function  of  the  stress 
in  the  elastic-plastic  layer,  crfp 


\+v 

KP 


af-K-pd’ia 


EP 

Jiju'kk  ’ 


(A-5) 


where  v  is  Poisson’s  ratio  and  I<P  is  the  stiffness  of  the  elastic 
element  in  this  network. 

Assuming  isotropic  strain  hardening  (represented  by  the  spring, 
H,  parallel  to  the  slider  element),  the  yield  strength  is  a  function  of 
the  plastic  strain  and  temperature,  0, 

°y  =  (  (^iP)  pi’  ^) '  (A-6) 

In  this  expression,  (ef.p)pl  is  the  equivalent  plastic  strain  given  by 
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Appendix  A.  Strain  formulations 


where  d(e|p)pl  are  the  increments  of  plastic  strain  tensor. 

Using  von  Mises’  yield  criterion  [45],  the  rate-independent 
plastic  flow  is  given  by 


AA.  Swelling  strains 

The  swelling  strains  in  the  reinforced  membrane,  caused  by 
water  uptake,  are  calculated  using  an  empirical  relationship  based 
on  volumetric  swelling  [43], 


f(aij)  —  yn^ij^ij  al 


uy j 


(A-8) 


where  ay  is  the  yield  strength  of  the  material  and  Sy  are  the 
deviatoric  stress  tensor  components  defined  by 
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$ij  —  & ij  ^  (A- 9) 

Under  von  Mises’  yield  criterion  [45],  the  material  deforms 
elastically  for/Oj/)  <  0  and  yielding  occurs  when 

/(<Ty,^')>0.  (A-10) 


A.4.  Strain  in  the  first  elastic-viscous  ( EV j)  network 


The  total  strain  in  the  elastic-viscous  network  (EVi)  is  the  sum 
of  linear-elastic  strain  in  the  spring  /CEVl  and  the  viscous  strain  in 
the  dashpot  element  given  by 


(»rv (#')„• 


(A-ll) 


where  the  subscript  el  refers  to  the  elastic  portion  and  vi  refers  to 
the  viscous  portion  of  the  strain.  Assuming  isotropic  Hooke’s  law 
for  the  linear-elastic  constitutive  response  of  the  spring  in  this 
network,  the  generalized,  three-dimensional  elastic  strain  as 
a  function  of  the  stress  crjjVl ,  is  given  by 


1+y 

^EVt 


^EV: 


—  d-aEV' 
I<ev,  y  kk  ' 


(A-12) 


For  the  dashpot  element,  the  Norton-Hoff  stress-strain  rate  law 
[46]  is  used,  the  generalized  three-dimensional  time-dependent 
form  of  which  is  written  as 


4  =  A(a f)",  (A-13) 

where  A  and  n  are  material  parameters. 
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